Strong field limit of black hole gravitational lensing 



O 

o 

(N 
D 



X 



V. Bozza*, S. Capozziello^, G. lovane-'-, G. Scarpetta^ 
Dipartimento di Fisica "E.R. Caianiello" , Universita di Salerno, Italy. 
Istituto Nazionale di Fisica Nucleare, Sezione di Napoli. 
(today) 

We give the formulation of the gravitational lensing theory in the strong field limit for 
a Schwarzschild black hole as a counterpart to the weak field approach. It is possible to 
expand the full black hole lens equation to work a simple analytical theory that describes 
at a high accuracy degree the physics in the strong field limit. In this way, we derive 
compact and reliable mathematical formulae for the position of additional critical curves, 
relativistic images and their magnification, arising in this limit. 
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I. INTRODUCTION 

The simplicity of the theory of gravitational lensing derives from some basic assumptions vifhich are 
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, satisfied by most physical situations y. Through the weak field and the thin lens approximations, the 

' whole theory was built and succeeded in explaining the phenomenology risen up to now. 

Yet, gravitational lensing must not be conceived as a weak field phenomenon, since high bending and 
looping of light rays in strong fields is one of the most well-known and amazing predictions of general 
^ ', relativity. The importance of gravitational lensing in strong fields is highlighted by the possibility of 

<sJ testing the full general relativity in a regime where the differences with non-standard theories would be 

^ \ manifest, helping the discrimination among the various theories of gravitation However, the com- 

W)' plexity of the full mathematical treatment and the difficulties of experimental or observational evidences 

are obstacles for these studies. 

Without the linear approximations, we have trascendent equations which are hard to be handled even 
. , numerically. 

' Several studies about light rays close to the Schwarzschild horizon has been lead: for example, Viergutz 

made a semi-analytical investigation about geodesies in Kerr geometry; in refs. 1^,^ the appearance 
of a black hole in front of a uniform background was studied. Recently, Virbhadra & Ellis faced 
the simplest strong field problem, represented by deflection in Schwarzschild space-time, by numerical 
techniques. The existence of an infinite set of relativistic images has been enlightened and the results 
have been applied to the black hole at the centre of the Galaxy. Later on, by an alternative formulation 
of the problem, Frittelli, Kling & Newman [Q attained an exact lens equation, giving integral expressions 
for its solutions, and compared their results to those by Virbhadra & Ellis. 

All these studies are affected by the high complexity in the investigation of the geodesies in the 
neighbourhood of the horizon and must resort to numerical methods to return valuable results. This 
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issue prevents any kind of general and systematic investigation of light bending in this region, since fair 
analytical formulae for the interesting quantities are absolutely missing. 

Starting from the black hole lens equation of ref. we perform a set of expansions exploiting the 
source-lens-observer geometry and the properties of highly deflected light rays. In this way, we manage 
to solve the lens equation and find analytical expressions for the infinite set of images formed by the 
black hole. This approach leads to extremely simple formulae which allow an immediate comprehension 
of the problem and a straightforward application to the physically interesting situations. This strong 
field approach can be surely considered as the direct counterpart of the weak field limit for its striking 
simplicity and reliability. On the other hand, this approach allows to investigate where the strong field 
limit effects become relevant. 

This paper is structured as follows: In Sect. 2, the black hole lens equation is presented and the first 
basic approximations are explained. Sect. 3 contains the calculation of the deflection angle. In Sect. 
4 the deflection angle is then plugged in the black hole lens equation to derive the relativistic images, 
their amplification and the critical curves. Sect. 5 contains a discussion of these results. Conclusions 
are drawn in Sect. 6. 



II. THE BLACK HOLE LENS EQUATION 

The geometrical configuration of gravitational lensing is shown in Fig. 0. The light emitted by the 
source S is deviated by the black hole L and reaches the observer O. Here, as a black hole, we mean 
any compact object having a radius comparable to its Schwarzschild radius, so that even very compact 
object which has not undergone a full gravitational collapse would work in the same way; P is the angular 
position of the source with respect to the optical axis OL and 9 is the angular position of the image seen 
by the observer. It is important to stress that the closest approach distance xq does not coincide with 
the impact parameter 6, unless in the limit of vanishing deflection angle a. 

By inspection of Fig. |l|, it is possible to write a relation among the source position, the image position 
and the deflection angle a. 

Dt q 

tan/3 = tan6l- — ^[tan6' + tan(a-6')]. (1) 
Dos 

This is what is called the full lens equation. Given a source position (3, the values of 6, solving this 

equation, give the position of the images observed by O. 

In the weak field limit, several standard approximations are performed. The tangents are expanded 

to the first order in the angles since they are, at most, of the order of arcsec. The weak field assumption 

4.GM 

reduces the deflection angle to — — . Then the lens equation can be solved exactly and two images are 

found: one on the same side of the source and one on the opposite. Their separations from the optical 
axis are of the order of the Einstein angle 



4GM Dls ^2) 
DosDoL 

It is easy to see that, in most relevant cases, these images are formed by light rays passing very far from 
the event horizon, justifying the weak field approximation [Q. 

Now we turn to the study of gravitational lensing in strong field. From now on, all lengths will 

2GM 

be expressed in units of the Schwarzschild radius — ^ — . The deflection angle a contains the physical 

information about the deflector and must be calculated through the integration of the geodesic of the 
light ray H. Its integral expression, as a function of xq, is 
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a (xo) = / , =dx — IT. (3) 

- -V(-)' 

The next section is completely devoted to the achievement of a manageable expression for this quantity 
as a function of 9. 

When the light ray trajectory gets closer to the event horizon, the deflection increases. At some 
impact parameter, a becomes higher than 27r, resulting in a complete loop of the light ray around the 
black hole. Decreasing further the impact parameter, the light ray winds several times before emerging. 
Finally, for b — (3\/3)/2, corresponding to xq = 3/2, the deflection angle diverges and the light ray is 
captured by the black hole. For each loop we add to the light ray geodesic, there is one particular value 
of the impact parameter such that the observer is reached by the light coming from the source. So there 
will be an infinite sequence of images on each side of the lens. 

We shall put our attention on situations where the source is almost perfectly aligned with the lens. In 
fact, this is the case where the relativistic images are most prominent. In this case, we are allowed to 
expand tan/3 and tan^ to the first order. Some more words are needed for the term tan {a — 9). Even if 
9 is small, a is not small in the situations of our interest. However, if a ray of light emitted by the source 
S is going to reach the observer after turning around the black hole, a must be very close to a multiple 
of 27r. Writing a = 2mr + Aa„, with n integer, we can perform the expansion tan {a — 9) ^ Aa„ — 9. 
The lens equation becomes 

f}^9- ^Aa„. (4) 
^os 

As defined by (^, a is a positive real number, corresponding to a clockwise winding in Fig. |]. We 
have inserted it in without giving a sign. Taking a positive /3, this equation describes only images on 
the same side of the source {9 > 0). To obtain the images on the opposite side, we can solve the same 
equation with the source placed in —f3. Taking the opposite of these solutions, we obtain the full set of 
the secondary images. 



III. THE DEFLECTION ANGLE 



The deflection angle can be evaluated exactly, but its expression does not allow the resolution of 
(^. Anyway, it is possible to make some simplifying but very general approximations which reduce the 
deflection angle to an expression easier to handle. By this strategy, we shall calculate this fundamental 
quantity. All approximations are essentially based on the proximity of the closest approach distance xq 
to its minimum value which is 3/2. 

The integral in Eq.(0) gives the following result 



4F (0o,A)G(xo) 



(5) 



where 



G(xo) 



Xq (^-3 + 3xo - + 2xo + x^j 
\ (3 - 2xo) (l~xo + ^-3 + 2xo + xl^ 



(6) 



and 



(1 - Asin^ 



-1/2 



(7) 
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is an elliptic integral of first kind. The parameters 0o and A themselves are functions of xq: 




y/-3 + 2X0 +^ 



and 



2(-3 + 2xo) 



.To - ^-3 + 2x0 + xl 



3 — Xq + + 2xo + Tq 

Fig. ^ shows a plot of the deflection angle. This expression is, of course, too complicated to use in 
Eq. (^. However, we can see that a diverges when xq —* 3/2. On the other hand, we are interested just 
into small closest approaches, since they correspond to the high deflection angles producing relativistic 
images. If we let xq — 3/2 + e, we can search for the leading order term in the divergence when e — > 0. 

After some non-trivial expansions, we find that the leading order of the deflection angle is logarithmic 
in e, that is 

a^-21og ^ ^ -TT. (10) 

In Fig. ^, we plot the ratio of the exact deflection angle and its leading order expansion. Of course, 
the two functions coincide in the limit xo — > 1-5. Using our expansion up to xo = 1.55, the error we 
commit is about 1%. For our purposes, this is largely sufficient as we shall see later. 

The next step is to convert the dependence on xq into a dependence on the impact parameter b and 
then on 9. The study of light ray dynamics in Schwarzschild metric gives the relation B 



(11) 

For values of a^o close to 3/2, this relation can be solved perturbatively (or, equivalently, one can expand 
the exact root of the third degree polynomial equation in xq to the first non-trivial order), giving the 
relation 



Remembering that b ~ 9Do l , we can substitute this relation into Eq. (|l^) to retrieve a as a function 
of 0, that is 

a^-\og(^DoL-^-^+A, (13) 

with 

(5 + 3V3) , , 

yl=-log-i -^-7r = 2.109. (14 

^ 1944 ^ ' 

Eq. ( p^ is the highly simplified expression for the deflection angle we were looking for. It allows 
reliable calculations, as we shall see in the next section, where we insert it in the lens equation. 
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IV. IMAGES, MAGNIFICATION, CRITICAL CURVES 



According to the considerations done at the end of Sect. 2, what enters the lens equation is not the 
fuU deflection angle but its deviation from a multiple of 27r. First of all, we have to find the values of 9 
(denoted by 0°) such that 

a [el) = 2n7r. (15) 
Inverting this equation, with a given by (|l3|), we have 

„o _ 3^3 + 2e'4-2"'^ 



2Do 



L 



(16) 
3 \/3 

As expected, if we let the number of loops around the black hole tend to infinity, we have 6*^^ — — — — . 

9^ are the starting points for all the successive calculations. The offsets Aq;„ can be found by expanding 
a to the first order in A6'„ = 9 — 9'^: 

^-n = -^^A9^. (17) 
We have all the ingredients to solve the lens equation (^), which becomes 

/3-KUA^«) + (^^)a^«. (18) 

We can observe that Dql, D^s, Dqs are all much greater than unity (remember that all distances are 
measured in Schwarzschild radii) and e'^"^^ = 3.25 x 10~^. This means that the last term in Eq.dl^) 
prevails on the A6'„ at the second place in the rhs. Neglecting this term, we finally get the position of 
the n*'' image as 



qO 



e 



{P^ 91) Dos 



DlsDol 



(19) 



We see that, when (3 equals 0°, there is no correction to the position of the n*'' image, that remains in 
9^-^ simply. In this particular case, the image position coincides with the source position. It is worthwhile 
to note that the second term in (|l9| ) is much smaller than the first one. For practical purposes, 0° are 
already a good approximation for the position of relativistic images. 

The relation (^) can help us to estimate the error in the determination of the position of the images. 
We do this in the least favorable case that is for the first image, which is the most external and thus the 
farthest from the divergence we started from. The relative error for 9\ is 

A^, fe---^a\Aa ^^^^ 



91 V DoL 91) a 

For the first image we have a ~ 27r. In the previous section, we have estimated the relative error on a to 
be about one percent for xq ~ 1.55. By using b = 9Dol and Eq. (p^), we can find the closest approach 
distance for the first image to be 1.545. In conclusion, this estimate for Aa/a can be reasonably used. 
The error on 9\ is 

— / ~ ~ 8 X 10"^ (21 

9\ 3\/3 + 2e^-2- a ^ ' 

For the other images, the errors can be found to be orders of magnitude smaller. This estimate 
supports our results which thus prove to be extremely reliable. 
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The critical curves are defined as the points where the Jacobian determinant of the lens equation 
vanishes, that is 

^^ = 0. (22) 

The term — is always positive (there are no radial critical curves). Then the (tangential) critical curves 

are obtained for (3 = 0. As we have already solved the lens equation for each it is sufficient to put 
/9 = in Eq.(p^) to get the angles of these relativistic curves, that is 

l-"- (23) 

A source perfectly aligned to the black hole produces an infinite series of concentric rings with these 
radii. 

The magnification of the images (|l^) is nothing else but the inverse of the modulus of the Jacobian 
determinant already used for the critical curves. In the region of our interest, i.e. for small source angles, 
it is easy to calculate all the needed quantities from the lens equation. We do this by approximating the 
images by the angles 0°: 



e.~^^ e^-^'-^Dos- ^^^^ 



The second term is much higher than 1 for the same reasons expressed before. Then the magnification 
of the n"* image is 

„ 1 € _A-2n. (3V3 + 2e^-^»-)i^os 

- I , , Ti r - ~ - "5 T^TTRi f; • y^^^i 



2f3Dl^DLs 



This expression gives a magnification decreasing very quickly and then the luminosity of the first image 
dominates all the others. The amplification diverges for (3 ^ 0, confirming the fact that the possible 
detection of relativistic images is maximal for a perfect alignment of the source with the lens. In our 
approximation, the amplification of the images on the opposite side of the source is just the same of 
those on the same side. 

Another interesting quantity is the total magnification of the relativistic images. As these would be 
seldom resolved as single images, it is likely to see them as one image with a total fiux equal to the sum 
of the partial contributions coming from each image. Then, we just have to sum up the series 

oc 

Mtot = 2^M„, (26) 

n=l 

which is a geometrical one, and then 



e- (3V3 (l + e^'^) + 2e-^) Dps 



(27) 



To take into account the finite extension of the source, one must integrate over its luminosity profile. 
As a simple example, integrating the -g dependence in Eqs. (|2^) and ( pTj) over a uniform disk with 
angular radius f3s, we get 



1 / 1 _ 2Stgn[ps-^] 




{Ps-Pf , 
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where we have indicated the disk with radius (3s centered on {3 by D(j3,l3s) and 



£;(</.o,A)= j {l-Xsm^^if^ d(t>, (29) 



is the elhptic integral of the second kind. 

Then, the magnification of an extended uniform source of angular radius j3s can be obtained by 
substituting the 1//3 dependence in Eq. (|25| ) and Eq. ( ^7|) with the right member of Eq. (^8|) . 

V. DISCUSSION 

The formulae just derived provide a complete characterization of the two infinite sets of relativistic 
images surrounding a black hole or, in general, any compact object acting as a lens whose size is 
comparable with its Schwarzschild radius. They can be employed in all kinds of phenomenological 
calculations to test their detectability. 

Some general considerations can be done on the features expected by a candidate lensing system. 
First, we need a very compact massive object, possibly a black hole, in order to gain access to the 
region of strong field immediately outside of the event horizon. Moreover, the matter surrounding this 
compact object should be transparent to the wavelength of the radiation emitted by the lensed source, 
otherwise the photons would be absorbed before the light rays complete their loops and would not reach 
the observer. 

We recall that the amplification of each of the weak field images can be expressed (with distances in 
Schwarzschild units) as: 



for the small Ps we are interested in. We can observe that the dependence on (3 is the same of the weak 
field images. Then the relative importance of the relativistic images to the weak field ones is constant 
for high alignments of the source to the lens. 

The ratio of this quantity and the total magnification of the relativistic images (|2^) is 



fJ'tot B \ Dos J 



(31) 



where _B is a numerical coefficient; 



e 



A 



(3^3 (1 + e2^) + 2e^) 



B = ^ ^ -{ = 0.017 (32) 

Then we must expect relativistic images to be always very faint with respect to the weak field images, 
since this ratio goes as an astronomical distance elevated to 3/2. This fact rules out microlensing as a 
method for detection of relativistic images. 

The separation between the two sets of relativistic images is of the order of the Schwarzschild diameter 
of the compact object. Depending on the specific situation, the angle corresponding to this length is 
generally very small. A very massive black hole would surely help to separate these images each other. 
The mass of the black hole at the center of our Galaxy §] is beheved to be about 2.8 x 10*^ Mq. Its 
Schwarzschild angle would be about 15 microsecs, which could become accessible by VLBI experiments 
[^,0. In this case, at least the two sets of relativistic images of some rear source would be distinguish- 
able. With their very low magnification, depending on the alignment degree, the relativistic images of 
any rear source would cast a very hard but not impossible challenge for observational astronomy. 

Possible better candidates could be black holes at the centers of other galaxies lensing some compact 
source (e.g. Quasars) on the background, or black holes at the centers of globular clusters. 
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VI. CONCLUSIONS 



The expansion we have presented in this paper can be considered as the strong field gravitational 
lensing limit, in opposition to the usual weak field limit. Its starting point is the critical null geodesic 
followed by a massless particle captured by the black hole. We have expanded the deflection angle in 
this limit and used its leading order in the full lens equation. We have solved this equation and found 
two infinite sets of relativistic images. We have given a complete analytical description of these images 
with simple formulae for their position, magnification and for the critical curves. 

The striking importance of these images lies on the fact that they could provide a profound test 
of general relativity in its full regime. In this situation, it would be surely possible to distinguish 
among relativistic theories of gravitation (e.g. Brans-Dicke, induced gravity, etc). For this reason, it is 
imperative to perform similar expansions within alternative pictures and evaluate the differences with 
the standard results we have found. 

The scheme we have used for this calculation encourages its application to more complicated problems 
of general relativity such as charged rotating black holes or other relevant astrophysical objects where 
strong field general relativity is involved. 
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FIG. 1. Geometrical configuration of gravitational Icnsing. Light rays emitted by the source S are deflected 
by the lens L and reach the observer O with an angle 6, instead of f3. The total deflection angle is a. xo is the 
closest approach distance and b is the impact parameter. Do l is the distance between the lens and the observer. 
Dls is the distance between the lens and the projection of the source on the optical axis OL. Dos = Dql + Dls- 
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FIG. 2. Deflection angle as a function of the closest approach distance. 



10 




11 



